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ABSTRACT 

Using homogeneous CCD photometric data from the 105-cm Kiso Schmidt telescope covering 
a 50' x 50' field, we study the mass functions (MFs) of nine open clusters. The ages and Galac- 
tocentric distances of the target clusters vary from 16 - 2000 Myr and 9-10.8 kpc, respectively. 
The values of MF slopes vary from -1.1 to -2.1. The classical value derived by Salpeter (1955) 
for the slope of the IMF is T — —1.35. The MFs in the outer regions of the clusters are found to 
be steeper than in the inner regions, indicating the presence of mass segregation in the clusters. 
The MF slopes (in the outer region as well as the whole cluster) undergo an exponential decay 
with the evolutionary parameter t (— age/ relaxation time). It seems that the evaporation of 
low-mass members from outer regions of the clusters is not significant at larger Galactocentric 
distances. It is concluded that the initial mass function (IMF) in the anticentre direction of 
the Galaxy might have been steeper than the IMF in the opposite direction. A comparison of 
the observed CMDs of the clusters with synthetic CMDs gives a photometric binary content of 
-40%. 

Subject headings: Galaxy: open clusters and associations- Technique: photometric- Stars: luminosity 
function, mass function- 



1. Introduction 

The distribution of stellar masses that form in 
a star- format ion event in a given volume of space 
is called the Initial Mass Function (IMF) and to- 
gether with the star formation efficiency, the IMF 
dictates the evolution and fate of star clusters. 
Present estimations of the observed IMF do not 
constrain the nature of the IMF (see e.g. Kroupa 
2007). The universality of the IMF of open clus- 
ters is still an open question because elementary 
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considerations suggest that the IMF ought to de- 
pend on star- forming conditions (Larson 1998). 
Therefore it is important to find systematic varia- 
tions of the IMF with different star- forming condi- 
tions. Identifying these variations would allow us 
to study early cosmological events (Kroupa 2002). 

Open clusters possess many favorable charac- 
teristics for IMF studies e.g., clusters contain an 
(almost) coeval set of stars at the same distance 
with the same metallicity; hence, difficulties like 
complex corrections for stellar birth rates, life 
times, etc. associated with determining the mass 
function (MF) from field stars are automatically 
avoided. The observed MF of a star cluster can 
in principle be determined from the observed lu- 
minosity function (LF) using theoretical stellar- 
evolutionary models. Since the MFs of interme- 
diate/old age open clusters may be affected with 
time due to stellar as well as dynamical evolution- 
ary effects, we can only estimate the present day 
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mass function. 

In recent years, luminosity and mass functions 
have been determined for a number of open clus- 
ters using homogeneous photoelectric or CCD 
data and reliable cluster membership criteria (cf. 
Piskunov 1976; Sagar et al. 1986, 1988; Scalo 
1986, 1998; Kjeldsen & Fransden 1991; Phelps 
& Janes 1993; Massey 1995; Durgapal & Pandey 
2001; Pandey et al. 2005, 2007 and references 
therein). Although the average slope of the 
MF does not seem to be very different from the 
Salpeter (1955) value, Pandey et al. (2001, 2005) 
found that the nature of the MF in open clusters 
does not remain the same over the entire region 
of the cluster and the slope of the MF steepens as 
radial distance from the cluster center increases. 

The nucleus and the corona (extended region 
of the star cluster) are two main regions in open 
clusters (Kholopov 1969). The nucleus of a cluster 
usually contains relatively bright i.e. massive (> 
3 Mq) stars along with low mass stars (see e.g. 
Brandl et al. 1999); whereas the corona, which 
contains a large number of faint and low mass (< 
1 Mq) stars, has important bearing on studies re- 
lated to the mass function. Extensive studies of 
the coronal regions of clusters have not been car- 
ried out so far mainly because of non-availability of 
photometric data in a large field around open star 
clusters. Therefore, considering the importance of 
low mass stars in the coronae of star clusters, we 
have generated a homogeneous wide-field photo- 
metric data-base around 9 open star clusters using 
a 2K x 2K CCD mounted on a Schmidt telescope 
(Kiso, Japan), covering a ~ 50' x 50' field (Sharma 
et al. 2006, hereafter Paper I). In this paper we 
study the mass function and photometric binary 
contents in those nine open clusters. 

2. Observational data 

Wide field broad-band CCD photometric ob- 
servations of the clusters were carried out during 
2001 November 19-25 using the 105 cm Schmidt 
telescope of the Kiso Observatory. A detailed de- 
scription of the observations and data reduction 
is given in our previous paper (Paper I). Briefly, 
the CCD camera used a SITe 2048x2048 pixel 2 
TK2048E chip having a pixel size of 24 x 24 /zm. 
At the Schmidt focus (f/3.1) each pixel of the CCD 
corresponds to l."5 and the entire chip covers a 



field of - 50' x 50' on sky. 

Initial processing of the data frames was done in 
the usual manner using IRAI0 and ESO-MIDA^l 
data reduction packages. Photometry of cleaned 
frames was carried out using DAOPHOT Software 
(Stetson 1987). The PSF was obtained for each 
frame using several uncontaminated stars. The 
FWHM of the star images varied between 3" and 
4" from night to night. The observations were cali- 
brated by observing standard stars in SA95 (Lan- 
dolt 1992) having brightnesses 12.2 < V < 15.6 
and color indices 0.45 < (B - V) < 1.51. Calibra- 
tion of the instrumental magnitudes to the stan- 
dard system was done using the procedure out- 
lined by Stetson (1992). The errors become sig- 
nificantly larger (>0.1 mag) for stars fainter than 
V = 20 mag; therefore, the measurements below 
this magnitude are not reliable and have not been 
used in the present study. 

3. Luminosity/ mass function 

With the help of color-magnitude diagrams 
(CMDs) we can derive the observed luminosity 
function (LF) of probable main sequence clus- 
ter members and then the mass function (MF) 
using theoretical evolutionary models, for which 
we adopt those of Bertelli et al. (1994). The 
factors which influence the determination of LFs 
from the observations are the accuracy of clus- 
ter parameters, data incompleteness and field-star 
contamination. The estimation of these factors 
and their treatment are described in the following 
subsections. 

3.1. Reddening, distance and age of the 
clusters 

The cluster parameters, listed in Table 1, are 
derived using the CMDs as discussed in paper 
I. The CMDs for stars lying within the clus- 
ter regions show a well-defined and broad main- 
sequence (MS). Barring Be 62, other clusters man- 
ifest a uniform reddening in the cluster region. 
Since the error in magnitude estimation for stars 
with V <18 mag is < 0.05 mag, we can conclude 
that the presence of probable binaries and field 
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stars should be the main cause for broadening of 
the MS in these clusters. In the case of Be 62 
variable reddening in the cluster region along with 
the presence of probable binaries and field stars, 
should be the cause of the broad MS. 

The extinction towards the clusters Be 62, NGC 
1960, NGC 2301 and NGC 2323 was estimated 
using the (U — B)/{B — V) two color diagram, 
whereas in the case of the other five clusters NGC 
1528, NGC 2287, NGC 2420, NGC 2437 and 
NGC 2548, the extinction was estimated using the 
V/(B - V) or V/(V - I) CMDs. The reddening 
in Be 62 varies from E(B — V) m i n = 0.70 mag to 
E(B-V) max — 1-00 mag. The distances and ages 
of the clusters were obtained by visual fitting of 
the theoretical isochrones by Bertelli et al. (1994) 
for Z = 0.02 to the blue envelope of the observed 
MS except in the case of NGC 2420, where we used 
isochrones for Z = 0.008, as Lee et al. (2002) have 
reported Z = 0.009 for this cluster. The accuracy 
of the distance estimates is ~ 10 per cent, while 
that of age determination is about 20 percent. The 
estimated values of E(B — V), distance and age of 
the target clusters (cf. Paper I) are given in Table 
1 and have been used in further analysis. 

3.2. Radial extent of the clusters and field 
region 

In paper I, we have studied the radial extent 
and structure of these clusters. The center of 
the cluster was estimated by convolving a Gaus- 
sian kernel with the stellar distribution and tak- 
ing the point of maximum density as the center. 
Projected radial stellar density in various concen- 
tric circles was obtained by dividing the number 
of stars in each annulus by its area. The extent 
of the cluster l r c i is defined as the point where 
the radial density becomes constant and merges 
with the field-star density. Within the uncertain- 
ties, the King model (King 1962) reproduces well 
the radial density profiles (RDPs) of the clusters 
studied in the present work. As an example the 
RDP of the cluster NGC 1960 along with the fit- 
ted King profile is shown in Figure 1. For this 
cluster, the core radius 'r c ', defined as the radial 
distance at which the value of radial density be- 
comes half of the central density (cf. paper I), 
and cluster extent l r c { come out to be 3.2 ±0.4 ar- 
cmin (1.2 pc) and 14 arcmin (5.4 pc), respectively. 
The structural parameters obtained by fitting the 



King- model surface-density profile to the observed 
radial density-profile of main-sequence stars hav- 
ing V < 18 mag arc taken from paper I and given 
in Table 2. 

It is well established that clusters have extended 
regions (coronae). Field-star contamination in- 
creases considerably in the coronal region of the 
cluster. The present observations have been made 
in a wide field (50' x 50'); the region outside the 
cluster extent (1.5 x r c /) has been used to estimate 
the field-star contamination in the cluster region. 

3.3. Probable members and data incom- 
pleteness 

To study the LF/MF, it is first necessary to 
remove field-star contamination from the sample 
of stars in the cluster region. In the absence of 
a proper-motion study, we used a statistical crite- 
rion to estimate the number of member stars in the 
cluster region. On the basis of a single passband 
alone it is difficult to establish that a particular 
star is in fact a member of the cluster. Therefore, 
two passbands, such as V and /, are required to 
identify the cluster members. We used V/ (V — I) 
CMD to estimate the membership as well as the 
LF of the cluster. The contamination due to field 
stars is greatly reduced by selecting a sample of 
stars which are located near the well-defined MS 
as described by Pandey et al. (2001, 2005). The 
same envelopes were used for the V/ (V — I) CMD 
of the field region to estimate the contamination 
in the cluster region due to field stars. After nor- 
malizing the area we can find the number of field 
stars which are considered to be present per unit 
area in each magnitude bin. As an example, selec- 
tion of the MS sample in the case of NGC 1960 is 
shown in Figure 2. 

The photometric data may be incomplete due 
to various reasons e.g. crowding of the stars, de- 
tection limit etc. The incompleteness correction 
is necessary if we want to analyze the LF/MF of 
the stars in the cluster. To determine the com- 
pleteness factor we used the ADDSTAR routine 
of DAOPHOT II. This method has been used by 
various authors (cf . Pandey et al. , 2005 and refer- 
ences therein). Briefly, the method consists of ran- 
domly adding artificial stars of known magnitude 
and position into the original frame. The frames 
are re-reduced using the same procedure used for 
the original frame. The ratio of the number of 
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stars recovered to those added in each magnitude 
interval gives the completeness factor, CF, as a 
function of magnitude. 

In practice we followed the procedure given 
by Sagar and Richtler (1991) and added artificial 
stars to both V and / images in such a way that 
they have similar geometrical locations but differ 
in I brightness according to mean (V — /) colors of 
the MS stars. The luminosity distribution of arti- 
ficial stars is chosen in such a way that more stars 
are inserted into the fainter magnitudes bins. In 
all about 15 per cent of the total stars are added 
so that the crowding characteristics of the origi- 
nal frame do not change significantly (cf. Sagar 
& Richtler, 1991). To have satisfactory statistics 
for the determination of CF, a number of indepen- 
dent sets of artificial stars are inserted into a given 
data frame (see for e.g. Table 3). The minimum 
value of the CF of the pair thus obtained is used 
to correct the data for incompleteness (cf. Sagar 
& Richtler, 1991). As an example, the CF along 
with relevant information for NGC 1960 is given 
in Table 3. As expected the incompleteness of the 
data increases with increasing magnitude and in- 
creasing stellar crowding. 

The number of probable cluster members in the 
two sub-regions of the cluster was obtained by sub- 
tracting the contribution of field stars (corrected 
for data incompleteness) from each magnitude bin 
of the contaminated sample of MS stars (also cor- 
rected for data incompleteness). The statistics in 
the case of one of the clusters, NGC 1960, is given 
in Table 4. 

4. Results 

4.1. Mass function 

The MF is often expressed by the power law, 
iV(logm) oc m r , where the slope of the MF is 
given as: 

r = d log N (log m)/d log m 

where iV(logm) is the number of stars per unit 
logarithmic mass interval. The classical value de- 
rived by Salpeter (1955) for the slope of the IMF is 
T = —1.35. The main-sequence LF, obtained with 
the help of CMDs for two sub-regions of the target 
clusters, has been converted into a MF using the 
theoretical model of Bertelli et al. (1994). The 



resultant MF data for the cluster NGC 1960 are 
given in Table 5. The MFs of the target clusters 
for two sub-regions as well as for the whole clus- 
ter region are shown in Figure 3. In the specified 
mass range the MF can be represented by a single 
power law. The value of MF slopes T, obtained 
by using the least-square solution in the specified 
mass range has also been given in Table 6. For 
intermediate age clusters the mass ranges are in 
general - 1 - 3M Q but for Be 62 and NGC 1960, 
which are the youngest clusters in the sample, the 
mass range is slightly higher. Despite large errors 
in r values, Table 6 indicates that in seven out 
of nine clusters the values of T are steeper in the 
outer region as compared to that in the inner re- 
gion. Barring the cluster NGC 1960, the difference 
in the values of T for inner and outer regions is less 
than 3 a (<~ 1 — 2 a) . The steeper values of T in 
the outer regions may be attributed to mass seg- 
regation. There is evidence of mass segregation in 
some Galactic and LMC clusters, with higher mass 
stars preferentially located towards the center of 
the cluster (see e.g. Fisher et al. 1998; Pandey et 
al. 1992, 2001, 2005 and references therein, Ku- 
mar et al. 2008). 

To evaluate the degree of mass segregation in 
clusters, we subdivided the samples into two mass 
groups as indicated in Figure 4, which shows the 
cumulative distribution of MS stars MS M function 
of radius in two different mass groups. In the case 
of six clusters (Be 62, NGC 1528, NGC 1960, NGC 
2323, NGC 2420 and NGC 2437) Figure 4 reveals 
the effect of mass segregation in the sense that 
relatively massive stars tend to lie near the clus- 
ter center. In the case of Be 62, NGC 1960, NGC 
2323 and NGC 2420 the Kolmogrove-Smirnov test 
confirms the above mentioned mass segregation at 
a confidence level better than 99 per cent, whereas 
in the case of NGC 1528 and NGC 2437 the con- 
fidence level is better than 90 per cent and 95 per 
cent, respectively. 

4.2. Dynamical state of the clusters 

Observations of mass segregation in several 
young clusters in the Galaxy (e.g. Moffat 1970; 
Hcrbst & Miller 1982; Larson 1982; Sagar et al. 
1988; Pandey, Mahara & Sagar 1992; Hillenbrand 
1997; Raboud & Mermilliod 1998) as well as in 
the Magellanic Clouds (Fischer et al. 1998 and 
references therein) suggest that mass segregation 
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may be the imprint of the star formation process 
itself. On the other hand, if clusters had a uniform 
spatial stellar-mass distribution at the time of for- 
mation, the spatial stellar-mass distribution would 
change with time as clusters evolve dynamically. 
Because of equipartition of energy the low-mass 
stars would attain high velocity and move away 
from the cluster center, consequently higher con- 
centration of high mass stars towards the center of 
the cluster could be observed (cf. Mathieu 1985; 
Mathieu & Latham 1986; McNamara & Sekiguchi 
1986). To decide whether mass segregation is pri- 
mordial or due to dynamical relaxation, we have 
to estimate the dynamical relaxation time, TV 
the time in which the individual stars exchange 
sufficient energy so that their velocity distribution 
approaches that of a Maxwellian equilibrium. The 
dynamical relaxation time is given by: 

rp _ 8.9x10 s N 1 / 2 ^ 2 _ 
E ~ fh 1 / 2 log{0AN) ' 

where N is the number of cluster stars, Rh is the 
radius containing half of the cluster mass and m is 
the average mass of cluster stars (Spitzer & Hart 
1971). 

We have estimated the relaxation time 'TV for 
all the target clusters to decide whether the mass 
segregation discussed above is primordial or due 
to dynamical relaxation. The total number of MS 
stars and the total mass of the MS stars in the 
given mass range (see Tabic 7) are obtained with 
the help of the MF. This mass should be consid- 
ered as a lower limit to the total mass of the clus- 
ter. For the half-mass radius, we used half of the 
cluster extent (R c i) obtained from the optical data 
(see. Table 2). The cumulative distribution of all 
the stars, shown by dotted curve in Figure 4, in- 
dicates that the 50 per cent of the cluster stars lie 
within ~ 0.48±0.10 R c i, therefore half of the clus- 
ter extent seems to be a reasonable approximation 
for the half-mass radius. The values of various pa- 
rameters as well as the resultant Te for the target 
clusters are given in Table 7. 

A comparison of cluster age with its dynamical 
relaxation time in the case of intermediate/old age 
clusters (age > 10 8 yr; see Table 7) indicates that 
the former is greater than the latter, leading to the 
conclusion that dynamical evolution could also be 
the reason for the observed mass segregation. In 



the case of the young clusters Be 62 and NGC 
1960, the dynamical relaxation time is compara- 
ble to the age of the clusters hence the observed 
segregation in these clusters could be because of 
both the imprint of the star formation process and 
dynamical relaxation. 

4.3. Synthetic CMDs 

During the last decade, synthetic CMDs have 
been used to study various properties of clusters, 
e.g., the MF and the influence of unresolved pho- 
tometric binaries on the LF, etc. (cf. Sandhu et 
al. 2003 and references therein). By comparing 
the synthetic integrated luminosity-function (ILF) 
and synthetic color-distribution with the corre- 
sponding observed distribution, these authors es- 
timated the photometric binary content in three 
intermediate-age open clusters. Following the pro- 
cedure of Sandhu et al. (2003) we calculated the 
ILF andA(V-J) = (V-/), - (V-I) M s for each 
star, where (V" — /)* is the observed color of a star 
and (V — I)ms is the corresponding color of the 
MS. The A(V — I) frequency distribution of stars, 
in the specified magnitude range (cf. Table 9) of 
the statistically cleaned CMD was compared with 
the A(V — I) frequency distribution of the syn- 
thetic CMDs (for details see Sandhu et al. 2003). 
Statistically cleaned CMDs were obtained using 
the following statistical procedure. For a star in 
the V, (V — I) CMD of the field region, the nearest 
star in the cluster's V, (V—I) CMD within V±0.25 
and (V — I) ± 0.13 of the field star was removed. 
While removing the stars from the cluster CMD, 
the number of stars in each luminosity bin was 
maintained as per the completeness corrected LF 
(cf. Sect. 3.3). 

Figure 5 shows the statistically cleaned CMD 
for the cluster NGC 1960. The CMDs of the clus- 
ters in the anticentre direction of the Galaxy are 
strongly affected by the background population 
of the Norma-Cygnus arm (see e.g. Pandey et 
al. 2006). The effect of the background popu- 
lation can even be seen in the statistically-cleaned 
CMDs. Therefore to avoid contamination due to 
background population we limit the MS popula- 
tion towards the fainter end and the fainter end 
limit is given in Table 9. 

Figure 6 shows the comparison of the observed 
ILF of the cluster NGC 1960 with the best-fit syn- 
thetic ILF for various percentages of photometric 
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binary content along with the obtained value of 
the MF slope. We assumed that the mass-ratio 
(mass of secondary / mass of primary) varies in the 
range of 0.75 — 1.0. The results for the target clus- 
ters are given in Table 8, which indicate that the 
value of the MF slope, T, for the observed MF is in 
agreement (within the errors) with the value of the 
MF slope obtained for the synthetic CMDs with- 
out binary content. Table 8 also indicates that the 
true/intrinsic value of T' becomes steeper if the 
photometric binary fraction is higher. This result 
is in agreement with that obtained by Sandhu et 
al. (2003) in the case of three open clusters. 

The comparison of observed and synthetic dis- 
tributions along with \ 2 values is shown as an ex- 
ample in Figures 7 and 8 for the clusters NGC 1960 
and NGC 1528, respectively. The results are given 
in Table 9 which indicate an average detectable 
photometric binary content of <~ 30% - 40% in the 
present sample. Mermilliod & Mayor (1989) found 
25% — 33% spectroscopic binaries in open clusters. 
Aparicio et al. (1990) and Durgapal et al. (2001) 
reported > 25% and 10% — 20% photometric un- 
resolved binaries in the clusters. In case of the 
Pleiades cluster, Bouvier et al. (1997) reported a 
binary (wider,visual) frequency of about 28 ± 4% 
for G and K dwarfs. Using the infrared speckle 
observations of the Hyades cluster, Patience et al. 
(1998) found that ~ 40% stars are binary. Ma- 
son et al. (1998) have estimated a binary fraction 
(both spectroscopic/unresolved as well as visual) 
of 75% in clusters/associations, whereas Jeffries 
et al. (2001), for the cluster NGC 2516, found a 
photometric binary fraction of 26 ± 5% for A to 
M-type systems with mass ratio between 0.6 and 
1. 

5. Discussions 

The mass functions of two clusters, NGC 1960 
and NGC 2323 are significantly steeper (> 3 a 
level) than the Salpeter value, whereas the MFs 
of the clusters NGC 1528, Be 62 and NGC 2437 
are found to be steeper but with lower significance 
level (2.1, 1.6 and 1.6 a level). The MF slopes for 
the youngest clusters Be 62 and NGC 1960 are 
based on a wide range of mass, i.e. 1.1-11.2 Mq 
and 1.0-6.8 Mq, respectively, while for other clus- 
ters the MF is derived for a relatively narrow mass 
range. The mass function slopes for outer regions 



are always steeper than the slopes for the inner re- 
gions. Mass function slopes of the youngest clus- 
ters, namely Be 62 and NGC 1960, are found to be 
-1.88±0.34 and -1.80±0.14, respectively, which 
are comparable to the slopes obtained for interme- 
diate/old clusters of the sample. The MF slopes 
for NGC 2287, NGC 2301, NGC 2420 and NGC 
2548 arc found to be comparable to the Salpeter 
value (-1.35). Bonatto & Bica (2005) have esti- 
mated r values of — 1.5±0.2 and — 1.3±0.2 for the 
clusters NGC 2287 and NGC 2548, respectively. 
For cluster NGC 1528, in the mass range 1.12- 
2.82 Mq, Francic (1989) has estimated a steeper 
mass-function slope (—2.78 ± 0.31). The present 
study also indicates a steeper mass function slope 
for NGC 1528. In the case of NGC 2323, for the 
mass range 0.40 - 3.90 M Q , Kalirai et al. (2003) 
have reported T = — 1.94 ± 0.15 which agrees well 
with the present value (-2.01 ± 0.17). For NGC 
1960, the present mass- function slope is steeper 
than the value given by Sanner et al. (2000) i.e. 
-1.23 ± 0.17 in the mass range 0.72-9.4 Mq. 

In order to investigate the relationship between 
relaxation time and cluster age with dynamical 
evolution, and to estimate the corresponding ef- 
fects on MFs we calculate for each cluster the evo- 
lutionary parameter, t, which is defined as the 
ratio of the cluster age to the relaxation time, 
t = age/Ts- Table 7 lists the estimated values of t 
for each cluster. Figure 9 shows T as a function of 
t. The slope of the MFs for six clusters, obtained 
from synthetic CMDs for 0% binary are shown by 
open circles. Data for two clusters (NGC 1907 
and NGC 1908) have been taken from Pandey et 
al. (2007) and shown by triangles. Although the 
errors in T values are large, Figure 9 clearly shows 
a systematic decreasing trend in T with t, partic- 
ularly in the outer regions of the clusters, indicat- 
ing a exponential decay of T with t. Bonatto & 
Bica (2005) and Maciejewski & Niedzieski (2007) 
have also concluded the same. However, the de- 
pendence of T in the central region of the cluster 
on t does not show the same trend. The decreas- 
ing trend may be represented by an exponential of 
the form: 

r = r + e a ' T 

with T = -1.9 ± 0.10, a = -17.6 ± 7.4 (corre- 
lation coefficient <~ 0.8, reduced x 2 = 0.04) and 
r = -2.11 ± 0.10, a = -18.6 ± 7.9 (correlation 
coefficient ~ 0.8, reduced \ 2 — 0.04) for the whole 
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cluster and outer regions, respectively, which in- 
dicate that the deceasing trend of T with r in the 
whole/outer cluster region is significant. Whereas 
for the inner region we estimate Tq = — 2.27±0.18, 
a = 0.11 ± 0.25 (correlation coefficient <~ 0.5, re- 
duced x 2 — 0.22) which indicates no correlation 
between T and t in the inner region. 

The parameter Ar = Ti nner — T outer , where 
Tinner and T oute r are the MF slopes for the inner 
and outer regions, respectively, can reveal infor- 
mation about mass segregation. Maciejewski & 
Niedzieski (2007) reported no correlation between 
Ar and r but Ar increases with age in the case 
of clusters older than ~ 100 Myr. In Figure 10 we 
plot Ar as a function of r and age, which does 
not reveal any relation between Ar and age of 
the cluster, however it seems that Ar decreases 
systematically with increase in r. The decrease 
in Ar with t can be interpreted as evaporation 
of low mass stars from the outer region. Figure 
10 also indicates a systematic variation of Ar as 
a function of Galactocentric distance, in the sense 
that Ar increases with increase in the Galactocen- 
tric distance, indicating that evaporation of low- 
mass members from outer region of the clusters is 
not significant at larger Galactocentric distances. 
Here we would like to point out that a larger sam- 
ple is needed to get a conclusive view about vari- 
ation of Ar with the t, age and Galactocentric 
distance. 

To study the dependence of the MF on the core 
radius r c , cluster extent r c ;, Galactocentric dis- 
tance Rg and age of the star cluster, we used 
the values derived in paper I. To convert the dis- 
tance to Galactocentric distance, the Galactocen- 
tric distance of the Sun is taken as 8.5 kpc (Allen 
2000). Figure 11 shows the dependence of clus- 
ter MF slope r on r c , r c i, Rg and cluster age. 
The Salpeter value for T (-1.35) is shown as a 
straight line. Figure 11 indicates that clusters 
having core radii greater than ~ 1 pc and clus- 
ter radii greater than ~ 4 pc have steeper MF 
than the Salpeter MF. The difference varies from 
1.6 - 2.1 a (Be 62, NGC 2437 and NGC 1528) to 
> 3 a (NGC 1960 and NGC 2323). The cluster 
located at Rg ~ 9.5 — 10 kpc also show a steeper 
MF at 1.6 a to 3.2 a level. As these clusters are 
situated in the anticentre direction of the Galaxy, 
it can be suggested that the IMF might have been 
steeper towards the anti-center direction as com- 



pared to other directions in the Galaxy. T does 
not show any trend with the ages of the clusters. 

6. Summary and conclusion 

In this paper we studied MFs of nine open star 
clusters located in the anti-center direction of the 
Galaxy, using wide field CCD photometric data 
taken from the Kiso Schmidt telescope. The values 
of MF slopes vary from -1.1 to -2.1. The main 
conclusions of the study are as follows. 

(1) The youngest clusters Be 62 and NGC 1960 
in the present sample have steeper MF 
slopes, -1.88 ± 0.34 and -1.80 ± 0.14 re- 
spectively than the Salpeter value. The 
observed MF of Be 62 can be assumed as 
the initial mass function (IMF), since the 
dynamical relaxation time is longer than the 
age of Be 62. 

(2) Three intermediate-age clusters (NGC 1528, 
NGC 2323 and NGC 2437) have steeper 
mass-function slopes, whereas four other 
intermediate-age clusters have MF slopes 
comparable to the Salpeter value. 

(3) Most of the clusters of the present sample 
show the effect of mass segregation. Mass 
segregation in the case of the young clus- 
ter Be 62 indicates that mass segregation 
could be due to the star formation process it- 
self, whereas in the case of intermediate/old 
clusters the mass segregation can also be ex- 
plained on the basis of dynamical evolution. 

(4) The MF slope of the outer region/ whole 
cluster region is seem to be related to the 
dynamical-evolution parameter r. The MF 
slopes (particularly in the outer region of the 
cluster) undergo an exponential decay with 

T. 

(5) There is evidence for initial mass-segregation 
within the young clusters and decrease in Ar 
with t is interpreted as evaporation of low 
mass stars from the outer regions of the clus- 
ters. 

(6) It is found that evaporation of low-mass 
members from outer regions of the clusters is 
not significant at larger Galactocentric dis- 
tances. 



7 



(7) The clusters having larger core/cluster radii 
have relatively steeper MF slopes. At larger 
Galactocentric distances the MFs of the 
clusters are found to be steeper. We do 
not find any correlation between MF and 
age of the clusters. 

(8) The present analysis of the synthetic CMDs 
reveals a detectable photometric binary con- 
tent of about 30% - 40% in the intermediate 
age clusters. 
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Table 1: Parameters of the target clusters taken from paper I. To determine the Galacto- 
centric distances Rq to the clusters, a value of 8.5 kpc (Allen 2000) has been assumed for 

the Galactocentric distance of the Sun. 

Cluster 012000 ^2000 E(B — V) Log age Distance Rq 







(h: 


m:s) 




(d:m:s) 


(mag) 


(y r ) 


(kpc) 


(kpc) 


Be 62 


01 


01 : 


15.5 


63 


56 


17 


0.70-1.00 


7.2 


2.32 


9.98 


NGC 1528 


04 


15 : 


24.2 


51 


15 


23 


0.26 


8.6 


1.09 


9.48 


NGC 1960 


05 


36 : 


20.8 


34 


08 


31 


0.22 


7.4 


1.33 


9.82 


NGC 2287 


06 


45 : 


58.7 


-20 


44 


09 


0.01 


8.4 


0.71 


8.96 


NGC 2301 


06 


51 : 


46.4 


00 


27 


30 


0.03 


8.2 


0.87 


9.25 


NGC 2323 


07 


02 : 


47.4 


-08 


20 


43 


0.20 


8.0 


0.95 


9.23 


NGC 2420 


07 


38 : 


24.8 


21 


34 


30 


0.04 


9.3 


2.48 


10.76 


NGC 2437 


07 


41 : 


58.1 


-14 


49 


28 


0.10 


8.4 


1.51 


9.51 


NGC 2548 


08 


13 : 


42.9 


-05 


46 


37 


0.03 


8.6 


0.77 


9.02 



Table 2: Structural parameters of the target open clusters taken from paper I. 









Optical 


2MASS 


Cluster 


1 


b 


Core radius 


Cluster extent 


Core radius 


Cluster extent 




(degree) 


(degree) 


arcmin (pc) 


arcmin (pc) 


arcmin (pc) 


arcmin (pc) 


Be 62 


123.98 


1.10 


2.2 (1.5) 


10 (7) 


2.5 (1.7) 


12 (8) 


NGC 1528 


152.06 


0.26 


8.3 (2.6) 


15 (5) 


18.5 (5.9) 


24 (8) 


NGC 1960 


174.52 


1.07 


3.2 (1.2) 


14 (5) 


3.8 (1.5) 


21 (8) 


NGC 2287 


231.02 


-10.44 


1.4 (0.3) 


12 (3) 


12.7 (2.6) 


16 (3) 


NGC 2301 


212.56 


0.28 


1.9 (0.5) 


9(2) 


4.5 (1.1) 


20 (5) 


NGC 2323 


221.67 


-1.33 


6.5 (1.8) 


17 (5) 


6.7 (1.9) 


22 (6) 


NGC 2420 


198.11 


19.63 


1.4 (1.0) 


10 (7) 


1.3 (0.9) 


9(7) 


NGC 2437 


231.86 


4.06 


6.8 (3.0) 


20 (9) 


9.6 (4.2) 


25 (11) 


NGC 2548 


227.87 


15.39 


1.5 (0.3) 


8(2) 


2.4 (0.5) 


8(2) 
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Table 3: Variation of the completeness factor with MS brightness in different radial 
(r, in arcmin) regions for the cluster NGC 1960. Nf is the number of frames generated; 
nc a ,ncn a ,nf a , nc r ,ncn r ,nf r and cft c ,cft cn , cftf are numbers of star added, recovered and 
corresponding completeness factors in inner, outer and field regions respectively. 



Inner region 
(r < 3'.2) 



Outer region 



Field region 
(17'. 5 < r < 22'. 5) 



Range in mag Nf 





nc a 


nc r 


cft c 


ncn a 


ncn r 




nfa 


nf r 


cftf 


V band: 






















13.6-14.6 


4 


19 


19 


1.00 


240 


237 


0.99 


277 


274 


0.99 


14.6-15.6 


7 


30 


30 


1.00 


480 


467 


0.97 


552 


542 


0.98 


15.6-16.6 


10 


53 


46 


0.87 


740 


723 


0.98 


892 


878 


0.98 


16.6-17.6 


13 


62 


58 


0.94 


1153 


1109 


0.96 


1268 


1243 


0.98 


17.6-18.6 


17 


83 


73 


0.88 


1674 


1549 


0.93 


1841 


1739 


0.94 


18.6-19.6 


21 


153 


118 


0.77 


2462 


2229 


0.91 


2662 


2475 


0.93 


I band: 






















13.6-14.6 


4 


16 


16 


1.00 


256 


255 


1.00 


293 


291 


0.99 


14.6-15.6 


7 


27 


26 


0.96 


528 


520 


0.98 


599 


592 


0.99 


15.6-16.6 


10 


54 


49 


0.91 


808 


789 


0.98 


920 


909 


0.99 


16.6-17.6 


13 


71 


63 


0.89 


1168 


1122 


0.96 


1368 


1332 


0.97 


17.6-18.6 


17 


88 


76 


0.86 


1743 


1626 


0.93 


1888 


1778 


0.94 


18.6-19.6 


21 


143 


110 


0.77 


2609 


2310 


0.89 


2802 


2565 


0.92 



Table 4: Luminosity function for the two sub-regions and the whole cluster region in the 
case of NGC 1960. n c , nf are the numbers of stars (corrected for data incompleteness) in 
the sub-regions of the cluster and expected in the field, respectively, and n p is the number 
of probable cluster members. 



Range Inner region Outer region Whole region 



V mag 


n c 


n f 


n p 


n c 


nf 


n p 


n p 


09.6-10.6 


2 





2 


5 


2 


3 


5 


10.6-11.6 


7 





7 


12 


1 


11 


18 


11.6-12.6 


10 


1 


9 


24 


11 


13 


22 


12.6-13.6 


11 


1 


10 


46 


23 


23 


33 


13.6-14.6 


16 


2 


14 


71 


45 


26 


40 


14.6-15.6 


19 


5 


14 


133 


89 


44 


58 


15.6-16.6 


17 


10 


7 


244 


190 


54 


61 



Table 5: Mass function of the cluster NGC 1960. The number of probable cluster members 
(N) was obtained after subtracting the expected contribution of field stars in each magnitude 
range, log <j> represents log (diV/dlog to). 



Range 


Mass 


Mean 


Inner 


region 


Outer 


region 


Whole 


region 


V (mag) 


(M e ) 


logM 


N 


log 4> 


N 


log (j) 


N 


log <f> 


9.6-10.6 


6.82-4.97 


0.7703 


2 


1.1625 


3 


1.3386 


5 


1.5604 


10.6-11.6 


4.97-3.45 


0.6241 


7 


1.6458 


11 


1.8421 


18 


2.0559 


11.6-12.6 


3.45-2.35 


0.4621 


9 


1.7303 


13 


1.8900 


22 


2.0983 


12.6-13.6 


2.35-1.74 


0.3102 


10 


1.8861 


23 


2.2478 


33 


2.4046 


13.6-14.6 


1.74-1.42 


0.1985 


14 


2.2004 


26 


2.4692 


40 


2.6563 


14.6-15.6 


1.42-1.18 


0.1136 


14 


2.2393 


44 


2.7367 


58 


2.8566 


15.6-16.6 


1.18-1.01 


0.0385 


7 


2.0085 


54 


2.8958 


61 


2.9487 
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Table 6: Mass- function slope T for two sub- regions and for the whole cluster region in the 
given mass range, a is the standard deviation of the slopes. 



Cluster 


Mass range 


Mass function slopes (T ± a) 






(M ) 


Inner region 


Outer region 


Whole cluster 


Be 62 


11.17-1.14 


-0.89 ±0.17 


-2.10 ±0.74 


-1.88 ±0.34 


NGC 


1528 


2.55-0.73 


-1.96 ±0.42 


-2.17 ±0.43 


-2.10 ±0.35 


NGC 


1960 


6.82-1.01 


-1.25 ±0.24 


-1.99 ±0.15 


-1.80 ±0.14 


NGC 


2287 


2.70-0.83 


-1.35 ±0.86 


-1.22 ±0.27 


-1.22 ±0.19 


NGC 


2301 


2.78-0.82 


-0.85 ±0.33 


-1.56 ±0.54 


-1.34 ±0.32 


NGC 


2323 


4.22-0.67 


-1.69 ±0.09 


-2.28 ±0.31 


-2.01 ±0.17 


NGC 


2420 


1.44-0.67 


-0.93 ±0.32 


-1.50 ±0.56 


-1.30 ±0.39 


NGC 


2437 


3.51-1.02 


-1.72 ±0.13 


-2.30 ±0.62 


-2.03 ±0.42 


NGC 


2548 


2.46-0.82 


-1.11 ±0.85 


-1.02 ±0.36 


-1.12 ±0.70 



Table 7: Various parameters for the target clusters in the given mass range (see Table 6) 
used for calculating the dynamical-evolution time Te- 



Cluster 


Number of 


Half mass 


Total mass 


Average mass 


Age Dynamical time 


T 




stars 


radius Rh (pc) 


(Me) 


m (Mq) 


(Myr) T E (Myr) (age/T E ) 


Be 62 


143 


3.4 


312 


2.18 


16 26 


0.6 


NGC 1528 


465 


2.4 


530 


1.14 


400 29 


13.8 


NGC 1960 


232 


2.7 


431 


1.86 


25 22 


1.1 


NGC 2287 


102 


1.2 


142 


1.39 


250 6 


41.7 


NGC 2301 


120 


1.1 


165 


1.38 


160 6 


26.6 


NGC 2323 


695 


2.3 


803 


1.16 


100 31 


3.2 


NGC 2420 


752 


3.6 


607 


0.81 


2000 75 


26.7 


NGC 2437 


647 


4.4 


1026 


1.59 


250 69 


3.6 


NGC 2548 


45 


0.9 


61 


1.36 


400 4 


100.0 


Table 8: Slope of the mass function Y obtained from synthetic CMDs for various percentages 




of binary content. 












Binary % 


NGC 1528 


NGC 1960 


NGC 2287 


NGC 2301 


NGC 2323 NGC 2420 







-2.02 ±0.13 


-1.85 ±0.12 


-1.31 ±0.21 


-1.32 ±0.16 


-2.10 ±0.13 -1.06 ±0.19 




10 


-2.14 ±0.20 


-1.79 ±0.16 


-1.36 ±0.18 


-1.57 ±0.17 


-2.16 ±0.10 -1.70 ±0.27 




20 


-2.43 ±0.17 


-2.01 ±0.11 


-1.51 ±0.18 


-1.65 ±0.15 


-2.27±0.17 -1.92 ±0.41 




30 


-2.51 ±0.23 


-2.28 ±0.18 


-1.52 ±0.15 


-1.72 ±0.11 


-2.37 ±0.13 -2.20 ±0.45 




40 


-2.67 ±0.30 


-2.29 ±0.20 


-1.62 ±0.28 


-1.88 ±0.25 


-2.53 ±0.13 -2.48 ±0.68 




50 


-2.68 ±0.26 


-2.29 ±0.20 


-1.70 ±0.16 


-1.90 ±0.19 


-2.66 ±0.22 -2.58 ±0.59 




Table 9: Binary fraction 


in various clusters. The expected error in estimation of binary 




content is <~ 


10 per cent. 












C luster 


V range Mass-range Photometric binary 










(mag) 


(M e ) 


content (%) 








NGC 1528 


11-17 


0.8 -3.0 


40 








NGC 1960 


10-17 


0.9 -6.1 


30 








NGC 2287 


09-16 


0.8 -3.2 


30 








NGC 2301 


10-17 


0.7 -3.1 


40 








NGC 2323 


10-17 


0.8 -3.9 


75 








NGC 2420 


15-19 


0.7 -1.3 


40 









12 



6 



5 - 




1 1 1 

5 10 15 20 25 

Radius (arcmin) 



Fig. 1. — The variation of stellar surface density for stars having V < 18, as a function of radial distance in 
NGC 1960. Continuous curve shows a least square fit of the King (1962) profile to the observed data points. 
The error bars represents ±y/N errors. The dashed line indicates the density of field stars. 




0.8 1.6 2.4 0.8 1.6 2.4 



(V-I) 

Fig. 2. — V/(V — I) diagrams for stars in the NGC 1960 cluster and field region. The slanted lines envelop 
the probable main-sequence stars. 
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Fig. 3. — A plot of the mass functions for two sub-regions and the whole-cluster region of the clusters Be 
62, NGC 1528 and NGC 1960. Log <p represents log (dTV/dlog m). The error bars represent ±y/N errors. 
Continuous curves show a least square fit for the given mass range. 



Inner region 




Outer region 



Whole region 



0.1 0.2 0.3 






i 


NGC 2301 

i 





• 1 ' 1 1 
















. NGC 2323 
. i.i.i 


i 







-0.2 -0.1 0.0 



0.1 0.2 0.3 0.4 0.5 0.6 

Log m ( M & ) 







NGC 2301 





















NGC 2323 









-0.2 -0.1 0.(1 0.1 0.2 0.3 0.4 (1.5 (1.6 



Fig. 3.— Continued - Same as before but for the clusters NGC 2287, NGC 2301 and NGC 2323 
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15 




Fig. 4. — The cumulative radial distribution of stars for different mass intervals. Dotted curves show the 
cumulative radial distribution of all the cluster stars. 



16 



14 



16 



18 



NGC 1960 




1.5 



Fig. 5. — Statistically cleaned CMD for the cluster NGC 1960. The solid curve represents the isochrone for 
14 Myr for Z=0.02 by Bertelli et al. (1994) adjusted for the distance and reddening of the cluster. 



300.0 



200.0 



100.0 



200.0 - 



100.0 



! binary=0% i/ '. 

- r=-i.85 / 


i i i i 1 i i i i 1 i i i i 1 i ii i 
: binary=10% ]/ '. 

- r=-i.79 / - 
yf\\ i i i i i i i i i,i i 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ii i 

: binary=20% 1/ ; 

- r=-2.oi / 
ktTT| i i i i i i i i i i 


ktTi | i i i i i i i ih 
'. binary=30% V - 

- r=-2.28 y 


binary=40% / - 
: T=-2.29 / ; 


i binary=50% }/ - 
- T=-2.29 7 - 



12.0 14.0 16.0 



12.0 



14.0 
V 



16.0 



12.0 14.0 16.0 



Fig. 6. — Synthetic integrated luminosity functions (ILFs) (continuous curves) of the cluster NGC 1960 are 
compared with the observed ILFs (filled circles). The value of binary content assumed and the values of 
mass function slopes obtained for the best fit are also given in the figure. The error bars represent ±>/~N 
errors in the observations. 
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-0.4 -0.2 0.0 0.2 0.4-0.4 -0.2 0.0 0.2 0.4-0.4 -0.2 0.0 0.2 0.4 
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Fig. 7. — Observed (filled points with error bars) and synthetic A(V-I) distributions (continuous curve) 
NGC 1960. 
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Fig. 9. — Variation of T as a function of r (age/TE) for the whole, outer and inner regions of the clusters. 
Filled and open triangles represent the slopes of the MFs for the clusters NGC 1907 and NGC 1912 taken 
from Pandey et al. (2007). The slope of the MFs for six clusters obtained from synthetic CMDs for 0% 
binary are shown by open circles. Continuous curves show a least square fit to the exponential law (see text). 
The data point for NGC 1912 (open triangle) is not included in the fit. 
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Fig. 10. — Variation of Ar {Tinner — r oute r) as a function of r (age/Te), log (age) and Galactocentric 
distance (Rg)- 



20 



-2 



1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 


| 1 1 1 | 1 1 1 | 1 1 1 | 1 1 1 

^ 1 ' 




1 — 1 

1 • 1 

1 • 1 

, , , 1 , , 




i , , , i , , , i , , , i , , , 



-1 



-2 



-3 



1 2 

r c (pc) 



4 6 

r cl (P°) 



9.5 10 



10.5 



7.5 



8 8.5 
log (age) 



10 



I 


i i i 


1 i 1 i 1 i 

, U- 


— 

I 1 1 1 1 


i 


i , i , i 


l 1 1 1 

1,1,1 


1 



9.5 



Fig. 11. — Dependence of the MF slopes (whole cluster) on core radius r c , cluster extent r c i, Galactocentric 
distance Rq and age of the clusters. The horizontal lines represents the value of Salpeter slope (r = —1.35). 
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